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We discuss how to extract renormahzed from bare Polyakov loops in SU (N) lattice gauge theories 
at nonzero temperature. Single loops in an irreducible representation are multiplicatively renormal- 
ized, without mixing, through mass renormalization. The values of renormalized loops in the four 
lowest representations of SU{3) were measured numerically on small, coarse lattices. We find that 
in magnitude, condensates for the sextet and octet loops are approximately the square of the triplet 
loop. This agrees with a large A'^ expansion, where factorization implies that the expectation values 
of loops in adjoint and higher representations are powers of fundamental and anti-fundamental loops. 
The corrections to the large relations at three colors are greatest for the sextet loop, ~ 1/A'^, and 
are found to be < 25%. The values of the renormalized triplet loop can be described by a matrix 
model, with an effective action dominated by the triplet loop: the deconfining phase transition for 
A'^ = 3 is close to the Gross-Witten point at A'' = oo. 



PACS numbers: 

I. INTRODUCTION 

In a SU{N) gauge theory, 't Hooft showed that the ex- 
act order parameter for deconfinement is a global Z{N) 
spin [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 
17, 18, 19, 20, 21, 22, 23]. The global Z{N) symme- 
try arises topologically [12] from the center of the gauge 
group, and is exact in a pure gauge theory. Quarks carry 
Z{N) charge, and so break the gluonic Z{N) symmetry. 
Nevertheless, numerical results from the lattice, termed 
flavor independence [24], suggest that the gluonic Z{3) 
symmetry may be an approximate symmetry of QCD. 

At a nonzero temperature T, a gluonic Z{N) spin is 
constructed by starting with a thermal Wilson line, which 
wraps all of the way around in imaginary time. The trace 
of the thermal Wilson line is the Polyakov loop [2] , 

^JV-^trLw, (1) 

and is gauge invariant. This is the trace of the propa- 
gator for an infinitely massive, test quark; the subscripts 
denote that the test quark is in the fundamental repre- 
sentation, of dimension N . 



As a gauge theory is heated, deconfinement occurs 
above a temperature T^. The confined phase is Z{N) 
symmetric, so the expectation value of the fundamental 
loop, which has unit Z{N) charge, vanishes below T^. 
The gluon spin condenses in the deconfined phase, above 

(M =e** KMI y^O , e*'^^ = l , T>Td, (2) 

and thereby breaks the global Z{N) symmetry sponta- 
neously [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 
17, 18, 19, 20, 21, 22, 23]. 

To compute near T^, it is necessary to employ numer- 
ical simulations on the lattice [24, 25, 26, 27, 28, 29, 
30]. The difficulty is that the expectation value of the 
Polyakov loop is a bare quantity, and so suffers ultraviolet 
divergences. This is due to an additive mass shift which 
the test quark undergoes with a lattice regularization. In 
four spacetime dimensions, the mass divergence for a test 
quark is linear in the ultraviolet cutoff, proportional to 
the inverse of the lattice spacing, a. This mass renormal- 
ization affects the expectation value of the bare Polyakov 
loop as the exponential of a divergent mass, mf}'", times 
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the length of the path [31, 32, 33, 34, 35, 36, 37]: 
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(3) 



For a Polyakov loop, the length of the path is 1/T. 

Gervais and Neveu [31], Polyakov [32], and others 
[33, 34, 35, 36, 37] established that Wilson lines are renor- 
malizable operators. Wc review these results in sec. II, 
applying them to a single Polyakov loop in an ar- 
bitrary, irreducible representation, TZ. A renormalized 
loop, is formed by dividing the bare loop by the ap- 
propriate renormalization constant, Zn: 



Zn 



Ztz = exp 



T 



(4) 



This is a standard type of mass renormalization [38] ; for 
example, in perturbation theory am^^ is a power series 
in the coupling constant. The only unusual feature is 
that because the Wilson line is a nonlocal operator, the 
renormalization constant depends upon the length of the 
path; in general, the renormalization constant for a path 
of length £ is = exp(— rn^"£). 

The real problem is how to extract the divergent 
masses non-perturbatively. In this paper we suggest a 
way of doing this. Consider a set of lattices, all at the 
same physical temperature, T, but with different values 
of the lattice spacing, a. Since the number of time steps, 
Nt — l/(aT), changes between these lattices, the diver- 
gent mass, arrij^" , follows by comparing the values of the 
bare Polyakov loops. This assumes that, as in perturba- 
tion theory, am*" is a function only of the temperature, 
and not of the lattice spacing, as a ^ 0. Given the renor- 
malization constant Z-ji, the renormalized loop ijz then 
follows from (4), up to corrections at finite lattice spacing 
~ aT. 

In an asymptotically free theory, at high temperature 
the vacuum is trivial in perturbation theory, as the ther- 
mal Wilson line is a Z{N) phase times the unit matrix, 
Ljv ^ c^'^Iat (2). After suitable normalization, the ex- 
pectation value of any renormalized loop approaches one 
at high temperature. 
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In sec. Ill wc present measurements of bare and renor- 
malized loops obtained through numerical simulations 
in a pure SU{'i) lattice gauge theory. Polyakov loops 
in the four lowest representations were measured, al- 
though we only found significant signals for three: the 
fundamental, the symmetric two-index tensor, and the 
adjoint representations. For three colors, these are the 
triplet, the sextet, and the octet representations, respec- 
tively. These loops were measured at temperatures from 
« .5Trf 3Trf. Numerically, we find that in all repre- 
sentations, the divergent masses am!^ are positive, so 
the bare loops vanish in the continuum limit, Nf oo. 



The values of renormalized loops appear to have a well- 
defined continuum limit and approximately satisfy the 
asymptotic condition of (5). 

An alternate procedure for computing the renormal- 
ized Polyakov loop was developed by Kaczmarek, Karsch, 
Petreczky, and Zantow [39]. They obtain Z-jz from the 
two point function of Polyakov loops at short distances. 
Their numerical values for the triplet Polyakov loop agree 
with ours at Td, but differ at higher temperatures; they 
did not consider higher representations. 

The most basic thing to consider is the size of the 
renormalized triplet loop. As we shall see, loops in higher 
representations are approximately given as powers of the 
triplet loop. We define a perturbative regime when the 
expectation value of the renormalized triplet loop is near 
one in magnitude. If the triplet loop is nonzero, but 
not close to one, then we have a deconfined, but non- 
perturbative, regime. At the transition, T = , both 
we and Kaczmarek et al. [39] find that the renormal- 
ized triplet loop, |(^3)|, is « .4; by ~ ^Td, we find it is 
w .9, while [39] finds ~ 1.0. This suggests that a pure 
gauge theory is, in some sense, perturbative from tem- 
peratures of « iTd on up. but not from Td to w iTd- 
This is in qualitative agreement with different resumma- 
tions of perturbation theory [7], all of which work down 
only to a temperature which is several times Td- It is also 
suggested by some features of the RHIC data [8]. 

Turning to loops in higher representations, we find 
that the expectation values of all renormalized loops are 
very small in the confined phase. Renormalized loops are 
nonzero above Td, with an ordering of expectation values 
as triplet, octet, and then sextet loop. Even though the 
octet and sextet condensates are smaller than that for 
the triplet, they are still significant. 

The apparently large values for the condensates of the 
octet and sextet loops are illusory. This is because when 
the fundamental loop condenses, that alone induces ex- 
pectation values for all higher loops. It is to these induced 
values that we must compare. 

This is clear in the limit of an infinite number of col- 
ors [10]. Migdal and Makeenko observed that in SU{N) 
gauge theories, expectation values factorize at large N 
[40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52]. Fac- 
torization is the statement that disconnected diagrams, 
with the most traces, dominate at large N . 

At infinite N, factorization fixes the expectation value 
of any Polyakov loop to be equal to powers of those for 
the fundamental and anti-fundamental, ij^ = (In)*, loop 
[41]: 
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The integers p+ and are determined from the Young 
tableaux of the representation TZ, using the composite 
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representations of Gross and Taylor [52, 53]. At any N, 
the overall phase is fixed, trivially, by the Z{N) charge 
of TZ, e-ji = P+ — modulo N. What is not trivial is 
the magnitude of the loop: at large A^, the term with the 
most powers of the fundamental loop dominates, with 
power p = p+ + P-. 

Lattice simulations with two colors by Damgaard and 
others [54, 56] showed that the bare adjoint loop is an 
approximate order parameter for dcconfincmcnt. At in- 
finite N, by factorization any renormalized loop serves 
as an order parameter for deconfinement, independent of 
its Z{N) charge. For example, consider the adjoint loop, 
with = p- = 1, and the loop for the symmetric two- 
index tensor representation, p+ — 2 and p^ — 0. While 
the adjoint has no Z{N) charge, and the two- index tensor 
charge two, modulo N, in magnitude both expectation 
values are ~ \{(!-n)\^ at large A^. 

We tested these large A^ relationships numerically for 
three colors. For each loop, we define the difference be- 
tween the measured loop and its value in the large N 
limit. The expectation value of the sextet difference is 
defined to be 
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and that for the octet difference as 



(8) 



(9) 



Of course there is some ambiguity in defining the dif- 
ference loops. One advantage of the above definitions is 
that they automatically vanish at both low and high tem- 
perature. In the confined phase, T < Td, the difference 
loops (nearly) vanish because all loops are (essentially) 
zero; at very high temperature, the difference loops van- 
ish because all loops approach one as T — > oo. 

The expectation values of the difference loops show 
interesting behavior. They vanish below T^, and spike 
down above T^, with a maximum at a temperature > T^. 
The spike for the octet difference is smaller, narrower, 
and closer to Td than the spike for the sextet difference: 
\(Sls) \ < .2, with a maximum at « l.lTd, while |(«6)| < 
.25, with a maximum at « l-STd. 

The magnitude of these expectation values are in ac- 
cord with a large A^ expansion. Corrections to the sextet 
difference are larger, 6£q ^ 1/N, than those for the octet 
difference, d£s ~ l/N"^. Thus our measurements of the 
values of renormalized loops give us a numerical estimate 
of just how good a large A^ approximation is for three col- 
ors. Corrections to the sextet difference, of ~ l/N, are 
found to be < 25% when A^ = 3. 

Although factorization tells us how to reduce conden- 
sates for higher loops to powers of that for the fun- 
damental loop, it does not tell us how the conden- 
sate for the fundamental loop changes with tempera- 
ture. Given the mean field relations satisfied by loops 
in higher representations, in sec. IV we consider a 
mean field theory for the fundamental loop itself. We 



consider a matrix valued mean field theory, or ma- 
trix model; this arises in a wide variety of contexts 
[41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52], includ- 
ing previous [48, 49, 54] and recent [14, 15, 16] work on 
the deconfining transition. 

The most general effective action for a matrix model 
of the deconfining transition involves a sum over loops 
in all representations. For three colors, we find that the 
lattice data for the renormalized triplet loop is approxi- 
mately described by a model whose action includes only 
the triplet loop. Over the range of temperatures stud- 
ied, the coiipling constant for the triplet loop is nearly 
linear in temperature. While the overall values of the 
sextet and octet loops are approximately described by 
this mean field theory, the difference loops are not. We 
categorize more involved models which might. 

It is interesting to consider what matrix models might 
apply to the deconfining transition for more than three 
colors. Consider the simplest possible model, where the 
action includes just the fundamental loop. At infinite 
A^, the solution follows from that of Gross and Witten 
[14, 15, 16, 43, 47, 49, 51, 54]. Kogut, Snow, and Stone 
showed that in this model, the deconfining transition is of 
first order, with a latent heat ~ [14, 15, 47, 49, 51, 54]. 
As a first order transition, the fundamental loop jumps, 
from zero to precisely one half: 



N 



N = oo . 



(10) 



This was also stressed recently by Aharony, Marsano, 
Minwalla, Papadodimas, and Van Raamsdonk [14]: at 
Tfi, this theory only deconfines halfway. 

It is striking that this special value at large A^ is close 
to that found from numerical simulations for the renor- 
malized, triplet loop; both we and Kaczmarek et al. [39] 
find 1(4)1 « .4 at T = T+. In the A^ = 3 matrix model 
closest to the Gross-Witten point, this value is not ^, 
but « .485 ± .001 [47, 54, 55]. The value in the Polyakov 
loop model, which we discuss shortly, is « .55 [19]. 

The first order transition at the Gross-Witten point 
is atypical. In ordinary first order transitions, masses 
are nonzero on either side of the transition [38]. Even 
though the value of the fundamental loop jumps at Td. 
at the Gross-Witten point both the string tension, and a 
gauge invariant Debye mass, vanish. This is only possible 
because of a transition, which is of third order in the 
matrix model coupling constant, at infinite A^ [43]. The 
Gross-Witten point is specific to infinite A^: at finite 
N, but > 3, in the matrix model deconfinement is an 
ordinary first order transition, with the string tension 
and the Debye mass nonzero at Td. 

For three colors, lattice simulations find a relatively 
weak first order transition, accompanied by a large de- 
crease in both the string tension and the Debye mass 
near Td, each by about a factor of ten [25]. The custom- 
ary explanation for this is that, as in the Potts model, 
three colors is near the second order transition known to 
occur for two colors [56]. We suggest that the deconfin- 
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ing phase transition for three colors is also close to the 
Gross- Witten point of infinite N; exactly how close can 
be categorized in a matrix model [55]. 

Aharony et al. [14] showed how at large N the Hage- 
dorn temperature can be computed when space is a very 
small sphere. They find that the deconfining transition 
is of first order if the Hagedorn temperature is greater 
than Td- It is tempting to think that the spikes which 
we found for the sextet and octet difference loops may 
be related to the Hagedorn temperature. If so, for three 
colors the Hagedorn temperature is tens of percent above 
that for the deconfining transition. 

It would be valuable to know from numerical simula- 
tions if the deconfining transition for more than three 
colors is close to the Gross-Witten point as well, or if 
that is unique to three colors. 

An Appendix gives a formal discussion of improved 
Wilson lines; on the lattice, these are related to smeared, 
stout links [57]. This may be of use for measuring 
Polyakov loops in higher representations. 

Our work was motivated by the Polyakov loop model, 
which postulates a relationship between the Polyakov 
loop — which was presumed to exist as a renormalized 
quantity and the pressure [18, 19, 20, 21]. In the end, 
we have more than expected: not just a renormalized 
Polyakov loop, but a good approximation to its poten- 
tial, in the SU{^) matrix model. 



II. BARE AND RENORMALIZED POLYAKOV 
LOOPS 

A. Traces of Wilson lines in imaginary time 

At a temperature T, the thermal Wilson line at a spa- 
tial point X, running in time from to r, is 

L7j(x,t) = V exp [ig ^S(x,r') dr'^ ; (11) 



we take the representation TZ to be irreducible. The nota- 
tion is standard: V denotes path ordering, g is the gauge 
coupling constant, r and r' are variables for imaginary 
time, T, r' : ^ 1/7", vector potential in the 

time direction, and are the generators of SU {N) in 
7^. 

The Wilson line is a SU{N) phase factor for TZ, and 
so is a unitary matrix. 



(12) 



d-R, is the dimension of the representation, and 1^^ the 
unit matrix in that space. 

The thermal Wilson line is proportional to the propa- 
gator of a "test" quark in the representation TZ [21, 33]. 
A test quark is one whose mass is so large that if you put 
it at a given point in space, x, it just sits there. The only 
motion of a test quark is up in imaginary time. While the 
test quark doesn't move in space, it still interacts in color 



space: through the Aharanov-Bohm effect, it acquires a 
SU{N) phase. 

To see this, form a covariant derivative in imaginary 
time, and define a propagator, Qqz, as its inverse: 



d_ 
J? 

It is easy computing the propagator in one dimension: it 
is just a step function, 9{t), times the Wilson line. 



gn{x,T)=e{T) LK(f,r). 



(14) 



Alternately, consider the path integral representation 
for the propagator of a particle with mass m in a back- 
ground gauge field; schematically, 

i^x'^^ds^ , (15) 

where = dx^/ds, and s is the path length; an exact 
form is given in [31]. In the limit of m ^ oo, this path is 
a straight line up in imaginary time, and this propaga- 
tor is the thermal Wilson line. Classically, the partition 
function in (15) is ^ cxp(— m£), where m is the bare 
mass, and £ is the length of the path. 

Under a gauge transformation in TZ, i^Tz{x, r), the Wil- 
son line transforms as 

Ltz{x, t) flliix, t)Ltz{x, t)Qtz{x, 0) . (16) 

As bosons, the gauge fields are periodic in imaginary 
time, with period 1/T. For the time being, we also as- 
sume that all gauge transformations are periodic in time, 
0,Ti{x, 1/T) = flTi{x,0). We relax that assumption later, 
but only in a way which affects the global symmetry, and 
not the local symmetry. 

For periodic gauge transformations, we can form a 
quantity which is locally gauge invariant by wrapping the 
Wilson line around in imaginary time, and then taking 
its trace: 



tr Lnix,l/T) 



(17) 



For the time being we follow the custom of mathematics, 
which is to work with traces which are not normalized. 
The trace is greatest when the Wilson line is the identity, 
= dn- 

The propagation of a test quark, at x, forward in imag- 
inary time generates the Wilson line in the fundamental 
representation, Ln{x,1/T). A test anti-quark is a test 
quark moving backward in imaginary time, so it gives the 
conjugate Wilson line, lij^ix, 1/T) = li\^{x, 1/T). 

Let us consider how to combine more test quarks and 
anti-quarks. To be gauge invariant, the Wilson lines must 
wrap around completely in imaginary time, so we drop 
the dependence on imaginary time, 1 /T. We also assume 
that all test quarks and anti-quarks are put down at the 
same point in space, and so drop the dependence on x 
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as well. So a test quark gives Ln, and a test anti-quark, 

Tt 

The general classification of representations is, for ar- 
bitrary iV, rather involved [52, 53]. We thus start with 
the lowest representations for general N. We then discuss 
some simplifications for N = 3. Finally, we show how one 
can easily classify all representations in the large N limit 
[52] . We use a notation where TZ is generally denoted by 
its dimension, d-ji. 

1. Simple examples 

The adjoint Wilson line is a test meson, constructed 
from a test quark and anti-quark. To combine the funda- 
mental and anti-fundamental Wilson lines into something 
with adjoint indices, we sandwich Ljv and Ljy between 
two SU{N) generators, t^, 

L^2_i = tr ^Ljv t% I'jv ; (18) 

the adjoint indices a, 6 = 1 . . . {N'^ — 1). The trace of the 
adjoint Wilson line is 

tr Ljv2_i = [tr Ljvi^ - 1 . (19) 

This follows from an identity on the t"'s, or directly from 
group theory. The product of the fundamental and anti- 
fundamental representations is the sum of the adjoint and 
identity, so the coefficient of jtr Ljvj^ in (19) is one. To 
check the —1 in (19), consider the case when the Wilson 
line is the unit matrix, L^r = Ijv; then the adjoint trace 
is its dimension, = N"^ — 1. 

The next representations are tensors with two funda- 
mental indices [53]. These represent the propagation of 
two test quarks up in imaginary time. Two Wilson lines 
can be put together in either a symmetric (-I-), or an anti- 
symmetric (— ), way. The Wilson lines for the two-index 
representations are 

^(jv2±Ar)/2 = 2 (^J^"'^^ ^ L-^I^Jv) • (20) 

Here i,j,k,l = 1...N are indices for the fundamental 
representation. The traces of the representations with 
two Wilson lines are then 

tr L(jv2±jv)/2 = ^ ((tr Lat)^ ± tr L^^ • (21) 

Checking the coefficients for = In, the dimensions of 
the representations are as indicated. The first term is a 
product of two Wilson lines, each of which wrap around 
once in imaginary time. The second term is one Wilson 
line which wraps around twice in imaginary time. 

2. Three colors 

For three colors, the fundamental representation is a 

triplet, 3, the adjoint representation is an octet, 8, and 
the symmetric two-index representation is a sextet, 6. 



Special to three colors, the anti-symmetric two-index 
representation is the anti-triplet, 3. To sec this, diago- 
nalize the triplet Wilson line by a local gauge rotation. 
After diagonalization, each element of L3 is a phase; the 
product of all phases is one: 

(exp(iai) \ 

exp(iQ!2) I . 

exp(-i(ai + a2)) J 

(22) 

Then it is easy to check that Lj^y^ -jY) /2 from (21) = tr Lg 
when = 3. Notice that (21) actually gives the anti- 
triplet loop, which is a fault of our notation. 

We mention one other representation for three colors. 
Consider a test baryon, composed of three test quarks. 
The symmetric combination of three fundamental Wilson 
lines gives the decuplet representation when N = 3. Its 
trace is 

tr Lio = ^ ((tr Ls)^ -|- 3 tr L3 tr -|- 2 tr L^^ . (23) 

Since L3 is a SU{3) matrix, 

det La = 1 



= ^ (tr L3)' - ^ tr L3 tr L2 + 1 tr , (24) 

which is Mandelstam's constraint [59]. Using it, we find 
that the trace of the decuplet Wilson line is 

tr Lio = tr Ls tr -M . (25) 
3. Large N 

The usual classification of representations is given us- 
ing Young tableaux [52, 53]. This is not very convenient 
for the large N limit, though. The reason is elemen- 
tary: Young tableaux involve the construction of tensors 
with fundamental indices. While of course complete, it 
does not naturally incorporate the symmetry between the 
fundamental and anti-fundamental representations. This 
is accomplished using the composite representations of 
Gross and Taylor [52]; we give an abbreviated summary 
which is sufficient for our purposes. 

Consider forming a state with test quarks and p_ 
test anti-quarks. This is done by combining funda- 
mental Wilson lines, and p- anti-fundamental Wilson 
lines, in all possible ways. To be gauge invariant, we 
must take traces of L^r and Ljy. We also have to remem- 
ber that Ljv is a unitary matrix, LjyLjv = Ijv, (12). This 
relation implies that all traces arc either traces of powers 
of Ljv , or traces of powers of Ljy , separately. Any mixed 
terms can be simplified using the unitary relation. 

The explicit construction of Wilson lines in different 
representations, as done above for the adjoint and the 
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two-index tensor, is unnecessary. Only gauge invariant 
quantities matter: these are traces of the Wilson line in 
different representations. For any TZ, by the Frobenius 
formula [52, 53] we can express the trace of a Wilson line 
in TZ as sums of products of traces of the fundamental 
Wilson line. 

We do this naively, by considering how to combine p+ 
Ljv's and p_ lJ^'s. The simplest thing is to take a trace 
of every Wilson line, separately: 

(tr L;v)^+ (tr L]^)p- . (26) 

In this term, each Wilson line wraps once, and only once, 
around in imaginary time. 

This is just the first term in a long series, though. To 
start with, we can consider a term where one Wilson line 
wraps around twice in imaginary time: 

(tr L;v)^+-' (tr L|,) (tr L]^)p- . (27) 

Continuing, a term where one Wilson line wraps around 
thrice: 

(tr L;v)^+-' (tr Li,) (tr L]^)^- , (28) 

and so on. We can do this as well for the anti- Wilson 
lines; a Wilson line which wraps around backward twice 
gives 

(tr L;v)^+ (tr Ljv)^'- tr {jJ^^f . (29) 

We continue in this fashion. The series continues down, 
generating fewer traces overall, as they are replaced by 
traces with higher powers of the Wilson and anti- Wilson 
line. In general, the operator tr represents a Wilson 
line which wraps around p+ times forward in imaginary 
time; tr (L^)^-, a Wilson line which wraps around p_ 
times backward. 

The series stops when we only have either one or no 
traces left. If p+ = p-, the scries stops at a constant, 
with no trace. If p+ ^ P-, the series stops at one trace, 

tr {iJ^Y- . (30) 

This term can obviously be reduced, as it is a Wilson 
line which goes forward p+ times, and then backward p- 
times. 

The above series of traces of Wilson lines represents 
the propagation of test quarks and anti-quarks. All we 
can do with test quarks is to take loops in imaginary 
time, so the only question left is how many times a test 
quark, or anti-quark, goes around in imaginary time. 

When p_ = these terms are the Schur functions [52, 
53]. When p_ ^ 0, these terms are the Schur functions 
for the composite representations of [52]. 

Consider now the powers of N at large N . We assume 
that in the deconfined phase, each trace gives a power 
of N: tr ~ AT for all q. The first term in (26) is of 
order ^ . Terms in (27) and (29) have one fewer trace 
overall, and so are ~ N^'^, down by one power of 1/7V 



relative to (26). Each time a Wilson line wraps aroimd 
an extra time in imaginary time, a possible trace is lost, 
which suppresses the term by ^ relative to (26). 

This assumes that the coefficients of all terms, (26-30), 
are of order one. For example, if the coefficient of (27) or 
(29) were ~ N relative to that of (26), then we couldn't 
conclude anything about the large N limit. Group the- 
ory tells us, however, that this is not so: for any N, all 
coefficients are numbers of order one [52, 53]. This is be- 
cause all representations are constructed by symmetriza- 
tion and anti-symmetrization operators [53]; the action 
of these operators depends on the number of indices, but 
not upon N. 

Consequently, at large A'' the trace of any representa- 
tion is dominated by the term where every Wilson line 
(or anti-line) wraps around only once in imaginary time. 
This is just because we maximize the number of possible 
traces, and so the powers of A^. At large A'' we denote 
representations asTZ = N /\p: 

tr LjvAp ~ (tr LnT^ (tr Lj^)f- , (31) 

p = p+ + P-. The notation is meant to be suggestive, as 
the dimension of this representation is ^ A^'' at large A''. 

The integers p+ and p- can be computed from the 
Young tableaux of the representation [52]. Denote the 
cohimns of the Young tableaux by the index /. and sep- 
arate them into two categories. If the number of rows in 
a column, r*, is < A''/2, then we leave the column alone, 
and refer to it as a fundamental column, with r!|_ = r* 
rows. If the number of rows in a column is > A^/2, 
then we turn it into a column with r!_ = N — r'^ anti- 
fundamental rows. Then p+ = S,;r^ is the number of 
boxes in all of the fundamental columns, and p- = Sjr^ 
is the number of boxes in all of the anti-fundamental 
columns. In the limit of infinite A^^, no other details of 
the Young tableaux matter; all that matters is the total 
number of boxes p+ and p- . 

In group theory, the distinction between fundamental 
and anti-fundamental columns appears awkward. It is 
essential to understand the large A^ limit, though. Con- 
sider, for example, a single anti-fundamental box in a 
Young tableaux: this is given by a column with N — 1 
rows. If we counted rows naively, this would give us A^ — 1 
powers of a trace. Clearly we should replace this by one 
anti-fundamental line, and one trace. Geometrically, if a 
Wilson line wraps forward around times in imaginary 
time, when > A^/2 it is better to replace it by a Wilson 
line going backward N — Vi times. 

While we do not need it now, we note that the Casimirs 
of the representations N A p were computed by Samuel 
[46] and Gross and Taylor [52]. While the dimension 
grows strongly, like A^^, the Casimir does not; at large 
A'', it is linear in N, proportional to p+ + p--. 

N 

Cnap = (p++P-) y +0(1) . (32) 

We shall see in sec (II. D) that this relation for the Casimir 
ensures that Polyakov loops, as computed to lowest order 
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in perturbation theory, satisfy factorization at large N, 
(58) and (55). 

B. Z(N) charges and confinement 

Previousfy, we required that gauge transformations be 
strictfy periodic in imaginary time, 't Hooft noted that 
this is not necessary. Consider gauge transformations 
which are periodic, in 1/T, onfy up to a constant. So 
as not to change the periodicity of the gauge fields, this 
constant gauge transformation must be equal to the iden- 
tity matrix times a phase, equal to an N*^ root of unity, 
(2). For the fundamental representation, aperiodic gauge 
transformations are 

Q.n{x, 1/T) = e'"^ 9.n{x, 0) , e'^^ = 1 . (33) 

This phase represents the center of the local SU{N) 
gauge group, which is a global group of Z{N): ^ = 
2nj/N,j = 0,l...{N-l). 

As the Z{N) phases commute with any element of 
the group, the gauge fields remain strictly periodic un- 
der such an aperiodic gauge transformation. The Wilson 
line, however, does change: for the fundamental repre- 
sentation, 

L;v ^ e'^ Ljv . (34) 

Wc define the Z{N) charge, or the A^-ality, of the fun- 
damental Wilson line to be one, bat = 1. The traces 
of Wilson lines in other representations transform under 
aperiodic gauge transformations as 

^ eienf . (35) 

the Z{N) charge e-jz is an integer. 

Due to the cyclic nature of Z{N), charge is only defined 
modulo N. If the fundamental has charge one, the anti- 
fundamental has charge minus one, which is equivalent 
to charge — 1. The simplest field with vanishing Z{N) 
charge is the adjoint. A baryon Wilson line, such as the 
symmetric combination of N fundamental indices, is also 
Z{N) neutral. At large A'^, the Wilson line in the N A p 
representation of (31) has Z{N) charge e-jz = p+ — p-, 
modulo A^. 

For three colors, the anti-triplet and sextet represen- 
tations have charge two, which is the same as minus one. 
As a test baryon, the decuplet Wilson line is Z{3) neu- 
tral. 

The confining phase, for T < T^, is characterized by an 
unbroken global Z{N) symmetry [1]. Hence the traces of 
Wilson lines with nonzero Z{N) charge vanish below T^, 

(tr Lk) = , T<Td , en¥=0- (36) 
Above Td, all traces develop expectation values, 

(tr L,^) ^ , T > Td , yen. (37) 



(Implicitly, we always assiimc that symmetry breaking 
occurs when a background Z{N) field is applied, and 
then allowed to vanish, in the appropriate infinite volume 
limit.) 

A priori, it is not obvious how the traces of Z{N) neu- 
tral Wilson lines behave in the confined phase, T < Td- 
Certainly, they must be nonzero at all temperatures. 
Even for T < Td, they are not protected by the Z{N) 
symmetry, and so are induced by quantum fluctuations 
at some level. For three colors, though, numerically we 
find that the expectation value of the renormalized octet 
loop is very small below Td, (63). This is natural in a 
matrix model, as discussed at the end of Sec. (IV. C). 

There is a counterpart to this in the behavior of large 
adjoint Wilson loops. At zero temperature, a fundamen- 
tal Wilson loop forms a string, with its expectation value 
the exponential of the string tension times the area. Ad- 
joint Wilson loops screen, so the adjoint string tension 
vanishes. Greensite and Halpern [40] show that at large 
N, the adjoint string breaks over distances which grow 
as ~ log(A'^). Similarly, the lattice finds that the adjoint 
string only breaks over large distances [29]. 



C. Prom traces of Wilson lines to Polyakov loops 

Traces of Wilson lines grow with the dimensionality of 
the representation. It is convenient to introduce a nor- 
malized quantity, which approaches one in the obvious 
perturbative limit. From the expectation value of the 
Wilson line in a given representation, we define the ex- 
pectation value of the Polyakov loop, £tz, as 

(in) = ^ (tr W) . (38) 
an 

The phase of the expectation of a given loop is fixed 
by the Z{N) symmetry. In the perturbative limit at high 
temperature, Ln e^^'^H-R, and 1-r e'^'^'^. 

Using normalized loops, instead of just traces, is most 
useful in considering the limit of a large number of colors. 
For example, the adjoint loop is 

1 N"^ ( I \ 

^^^-1 = N^^' ^^^-^ = A^ U^^'' ~ A^J • 

(39) 

The large A^ limit of this expression cannot be taken 
directly, since it involves an operator, j^Arj^, and a pure 
number. A large A^ limit can be taken by comparing ex- 
pectation values. This is especially easy at large A^, as 
then disconnected diagrams dominate, and all expecta- 
tion values factorize [40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 
50, 51, 52]. For instance, 

{M) = KMI' , ^ = 00 . (40) 

Using this, at large A^ the expectation value of the adjoint 
loop is, in magnitude, identically the square of that for 
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the fundamental loop, 

(€^._i)«|(MP+o(^) . (41) 

As it stands, this is a relationship between bare loops; 
we shall see, however, that it survives renormalization. 
For the two index tensor representations, 

2 

= 7^^(4±>L^) (42) 

Hence at iV = oc' the expectation values of the two-index 
tensor loops are the square of the fundamental loop, 

(V=±iv)/2>«(M' + 0(^1) . (43) 

Notice that the term at infinite N is the same for both 
the symmetric and the anti-symmetric representation; 
the difference only shows up in corrections in 1/A''. This 
generalizes to higher representations at large A'^. 

For the two index tensor representations, corrections 
in 1/N start with the expectation value of the operator 

In the deconfined phase, we consider the trace of any 
power of the Wilson line to be a number of order one, so 
the trace of is like that of Ljv, a number of order N. 
Overall, then, this operator is ^ 1/A'^. It is not surprising 
that corrections to the two-index loop are ~ 1/-/V, larger 
than the ^ l/N"^ for the corrections to the adjoint loop. 
These corrections for the two-index loop arise because 
the operator has Z(N) charge two, and so mix with the 
operator in (44). The adjoint loop is Z{N) neutral, and 
so can't mix with this operator. The analogous operator 
for the adjoint loop is tr LjyLAr, which by the unitary 
relation of (12) is a constant. 

The generalization to the N A p representations of the 
large N limit is immediate, 

which is (6) and (7). One advantage of using loops is that 
it is easy to check overall normalization: up to Z{N) 
phases, both sides approach unity in the perturbative 
limit. 

A systematic expansion in 1/A^ proceeds by including 
operators such as 



for integral q. This operator is ^ 1/N'^~^ in the decon- 
fined phase. As q grows, the number of such operators 
docs as well. We do not concern ourselves here with the 
development of a systematic expansion, and consider only 
the leading corrections in 1/A'. For three colors, our nu- 
merical simulations indicate that this might not be such 
a bad approximation. 

Even to leading order at large A^, some representations 
involve loops other than the fundamental. Consider a 
test baryon, composed of A^ fundamentals. Because of 
Mandelstam's constraint that Ljv is an SU{N) matrix, 
det Lat = 1 , the term with A" powers of the fundamental 
loop is part of the identity representation. Hence at large 
N, the loop for a test baryon behaves as 

{etest-baryon} = (^iv)^"' ^ ^ tr L^^ +0 (^1^ . (47) 

This is illustrated by the decuplet loop for three colors, 
(25). 

D. Renormalization of Polyakov loops 

With this lengthy introduction aside, we turn to the 
problem at hand, the renormalization of Polyakov loops. 

Remember how mass renormalization usually works, 
as for a scalar field, (/), in four spacetimc dimensions [38]. 
If the mass of the field is m, and its coupling X(j)'^, to one 
loop order the mass squared receives contributions 

a momentum cutoff A is used to regularize the integral. 

The structure at one loop order is generic to perturba- 
tion theory: there are two mass divergences, one propor- 
tional to a power of the cutoff, ~ A^, and the other, to 
a logarithm of the cutoff, ~ m^log(A/m). The power 
divergence is an additive shift in the bare mass, and for a 
scalar field is inconsequential: the parameters are tuned 
to be near a critical point, where the renormalized mass 
vanishes. On the other hand, the logarithmic divergence 
is physical, related to the anomalous dimension for the 
mass operator. A renormalization condition is required 
to fix the value of the renormalized mass at a given scale. 

Polyakov loops correspond to a test particle whose 
mass is taken to infinity, so their worldline is a straight 
line. This freezes out fluctuations in the time-like di- 
rection. This is obvious in perturbation theory: as 
/ Ao{x,T)dT always enters, only modes which are con- 
stant in r appear. Thus the mass divergence of a 
Polyakov loop in four spacetime dimensions is like that 
of a propagating particle in one less dimension, which 
is three. Similarly, the mass divergences of a scalar 
field in four spacetime dimensions, (48), are like those 
of Polyakov loops in five spacetime dimensions. 

The ultraviolet divergences of a Wilson line depend 
only upon the representation, and not upon (smooth) 
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details of the path. For the time being, let L^j denote 
any Wilson line in a representation TZ: we only assume 
that the path forms a loop, so that traces of are gauge 
invariant. The expectation value of the Wilson line has 
a mass divergence which depends upon the length of the 
loop, jC, as 

(tr Lk) ~ exp {-mtC) . (49) 

The exponentiation of mass divergences follows from 
the analysis of Gervais and Neveu [31]. Similar to (15), 
they rewrite the Wilson line as a propagator for a fermion 
which lives in one dimension, along the path of the loop. 
With lattice regularization, the additive mass shift which 
the Wilson loop undergoes, m^" , is no different from that 
which (non-gauged) propagating fields also experience, 
such as for massive quarks, or scalar fields, (f>. 

On the lattice, the exponentiation of mass divergences 
has been shown explicitly by Curci, Menotti, and Paffuti 

To develop insight into the divergent masses, we com- 
pute to one loop order. In four spacetime dimensions, 

mn ^ + Cng / "p" ~ + ■ (50) 

We have used a lattice, with lattice spacing a, to regu- 
larize the theory. The exact coefficient of 1/a in m*" 
depends on the details of the lattice discretization, but 
it is a positive, nonzero number of order one. In four 
dimensions, om*" is a power series in the coupling con- 
stant. 

For a straight Polyakov loop in four dimensions, this 
is the only divergence: there is no anomalous dimension 
for the corresponding mass. This is clear to any order in 
perturbation theory, and occurs because the mass diver- 
gence is like that of a particle which propagates in three, 
instead of four, dimensions. 

Loops can also have cusps [32, 33, 34, 35, 36]. In order 
to be periodic in imaginary time, the simplest example of 
a Polyakov loop with cusps has not one, but two, cusps. 
This is illustrated in Fig. 1, with cusps at r = and 
T = 1/(2T). These cusps reflect external probes which 
deflect the test particle at these points. As Polyakov 
loops, the expectation values of single loops with cusps 
only have nontrivial expectation values in the deconfined 
phase. 

A cusp generates a logarithmic singularity in four 
spacetime dimensions [32, 33, 34, 35, 36]. This is not 
proportional to the length, and so does not contribute 
to the divergent mass. A condition to fix the value of a 
renormalized loop with a cusp must be supplied, but this 
is standard. For example, in QCD loops with cusps are 
related to the Isgur-Wise function [36]. 

It is also interesting to consider loops in three, instead 
of four, spacetime dimensions [58]. The linear divergence 
of (50) is now logarithmic [31], 

mf " ~ + Cng' f'^ ^ ~ + Cng^ log ; (51) 



in three dimensions, the coupling constant has dimen- 
sions of mass. As the divergent mass depends logarithmi- 
cally upon the lattice spacing, a condition to fix the value 
of the renormalized loop must be supplied [60]. Loops 
with cusps do not have new ultraviolet divergences in 
three spacetime dimensions, although they do have power 
like infrared divergences. 

Defining wave function renormalization for a Wilson 
line of length C as 

= exp (- m^C) , (52) 
then the renormalized Wilson line is given by 

L-R. = ^ L-R. ; (53) 

as illustrated by the renormalization of Polyakov loops, 
(4). In the space of all SU{N) invariant tensors, the set 
of irreducible representations form a complete and or- 
thonormal basis [53] . As this basis is orthonormal, Wil- 
son lines in different reprc^sentations do not mix. Conse- 
quently, in different representations the divergent masses, 
m^", and so the renormalization constants Ztz, are in- 
dependent quantities. 

We discuss in the next section how to extract the di- 
vergent masses from lattice simulations using a straight 
Wilson line. In an Appendix we also discuss how the Wil- 
son line might be modified to alter the mass divergence. 

It is illuminating to compute the renormalized loops 
to one loop order. Then it is easiest using dimensional 
regularization, as then the divergent mass automatically 
vanishes. Following Gava and Jengo [3, 21], the leading 
correction arises after the Debye mass, m|j = 
is included by resummation. To lowest order, the correc- 
tion to the renormalized loop is 

(M-1 ^ -Cng' J f^2l ^ +Cng'^D (54) 
so that 

<^^'--^-°(»')- 

(In three spacetime dimensions, {i-jij — 1 ~ 
Cn{g'/T)\og{Tlg^) [60].) 

In four spacetime dimensions, the leading correction 
to the renormalized loop is positive. Thus in the limit of 
high temperature, the loop approaches one from above, 
and not from below. At first sight, this seems paradoxi- 
cal. Bare Polyakov loops are traces of SU [N) matrices, 
and so satisfy a strict inequality, \C.-]z\ < 1- For example, 
on the lattice this holds configuration by configuration. 
Instead, renormalized loops satisfy the renormalized con- 
straint, 

\^n\ < . (56) 
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Numerically, we find that the divergent masses are uni- 
formly positive, 

am^''>0, (57) 

for all representations, at all temperatures. If so, then 
the renormalization constant always vanishes in the 
continuum limit, and there is no constraint on the renor- 
malized loop. This condition is most natural: otherwise, 
Zjz diverges as a —> 0, so the renormalized loop must 
vanish. 

The renormalization of a constraint is also familiar 
from the non-linear sigma model in two spacetime dimen- 
sions [38]. If the sigma field is an SU{N) matrix, then 
like the bare Wilson line, the bare field is a unitary ma- 
trix. Because of wave- function renormalization, however, 
the renormalized sigma field satisfies a renormalized, and 
not a bare, constraint. 

In the large N limit, factorization holds. This implies 
that 

a " K{p++p-)a mj}" . (58) 

This is automatic to lowest order in perturbation the- 
ory, where am'^ ^ C-ji, remembering that the Casimirs 
satisfy (32). This also ensures that the perturbative ex- 
pression for the renormalized Polyakov loop, (55), is well 
behaved at large iV. 

McLerran and Svetitsky [4] used the expectation value 
of a loop to define the free energy of a test quark. If this 
test free energy is defined from the renormalized loop 
as J^Ti = — riog(|(^K)|), then while it is positive near 
Td, from (55) it is negative at high temperature, J^tz ^ 

-Wl0g(T)3/2. 

While the divergent masses depend upon the ultravio- 
let cutoff in a unremarkable manner, the renormalization 
constants are not like those of local operators. In four 
spacetime dimensions, the renormalization constants of 
local operators are independent of temperature. In con- 
trast, the renormalization constants of Polyakov loops 
are temperature cicpendcnt, but just because the length 
of the path for a Polyakov loop is 1/T. 

Renormalization implies that the only measurable 
quantities are single trac;es of Wilson lines. Consider the 

most general, gauge invariant combination of bare Wilson 

„+ 

lines possible. For example, start with tr L^^ , which rep- 

resents the propagator for a test quark, in the representa- 
tion TZf , q'^ times around some fixed loop in spacetime. 
Generically, we can take powers of this trace, and then 
multiply different powers of different traces together. We 
can also do the same with conjugate operators. By the 
character expansion [53], any such combination can be 
reduced to a linear sum over traces of single Wilson lines 
in different, irreducible representations: 

= ^CTC ^TC = . (59) 

n n 



Here all n's and g's are positive integers; the constants 
cqz, and the representations TZ which one must sum over, 
are determined by group theory [53]. Because irreducible 
representations form a complete basis over all SU{N) 
representations [53] , we can insist that only linear powers 
of Wilson lines appear on the right hand side. With 
a linear sum, renormalization is then just a matter of 
replacing bare by renormalized loops. 

Assuming that all m^" > 0, (57), so the Z-jz all vanish 
as a ^ 0, in the continuum limit only the identity rep- 
resentation survives. This is of no physical consequence: 
the physical quantities, the traces of renormalized Wil- 
son lines, are hidden in the corrections to this relation, 
which are exponential in 1 /a as a ^ 0. 

This was discovered numerically. To a high accuracy, 
we found that (l^sl^) « 1/9; corrections varied from ~ 7% 
for Ni = 4, to - .2% for iV^ = 10. This is because 
from (19), - 1/9 « 8Z8(4)/9; because of the octet 

renormalization constant, Zs, this is a small quantity. 

Previous work on the renormalization of loops at zero 
temperature concentrated on loops in the fundamental 
and adjoint representations, especially on the case of 
loops with cusps [35; 36]. The case of traces of lines 
which wrap around the same loop several times, or prod- 
ucts of such traces, was neglected. At nonzero tempera- 
ture, though, the natural loops to consider are those at 
the same point in space, wrapping around in imaginary 
time in all possible ways. As discussed, this is equivalent 
to the set of loops which wrap around in imaginary time 
just once, although in arbitrary representations. 

III. LATTICE MEASUREMENTS OF SU{3) 
POLYAKOV LOOPS 

A. General Method 

We turn to the case of three colors. Group theory 
tells us how bare loops are related, through expressions 
such as (39). After renormalization, we do not know 
how renormalized loops are related. Except at very high 
temperature, where we can use perturbation theory, (55), 
the only way to compute renormalized loops is through 
numerical simulations on the lattice. 

In this section we discuss how we extract renormalized 
Polyakov loops for the lowest, nontrivial representations 
of SU{3) color. Consider a lattice with Ng steps in each 
of the three spatial directions, and Nt steps in the time 
direction. At lattice spacing a, the physical temperature 
T = l/{aNt). As discussed in the Introduction, to ex- 
tract the mass divergence, we consider a series of lattices, 
all at the same temperature, but with different values of 
the lattice spacing, a, and so Nt] Ng/Nt is kept in fixed 
ratio. At a given value of T/T^, we assrune that the loga- 
rithm of the expectation value of a single, bare Polyakov 
loop can be written as a power series in l/Nt. 

- log {In) = f^^^Nt + + f^' ^ . (60) 
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In four spacctimc dimensions, f!j^ 



div 



am: 



div 

n 



(In three 



dimensions, the series is fSj^'" log(7Vt) + Z^-^" + f^^/Nf) 
Each of the Jtz'h is a power series in the couphng con- 
stant, g"^. On the lattice, this is a series in the bare 
couphng constant, and becomes, in the continuum hmit, 
a series in the renormahzed couphng constant. As such, 
the fns are functions only of the temperature divided by 
the renormalization mass scale; or equivalently, of T jT^. 
By comparing expectation values at the same tempera- 
ture, but different values of A^t, we can extract am^^. 
What remains is the renormalized loop in the continuum 
limit, 



£n) = exp(-/. 



ren\ 

n ) 



(61) 



There are also corrections at finite lattice spacing, /^*. 
Near the continuum limit, these effects begin at l/N^, 
with = Yl'^=i^j l^t~ ■ weak coupling, ci ^ g**; 
these are corrections to the one loop term on the lattice, 
after resumming the Debye mass. Corrections in C2 ^ 
presumably arise as lattice corrections to the continuum 
term, (61). 

As is common on the lattice, we work at a fixed ratio of 
Ng/Nf. Thus we implicitly assume that the dependence 
upon this ratio is negligible in the infinite volume limit. 
This can be studied analytically, but requires a careful 
treatment of the constant modes. A perturbative study 
is given by Heller and Karsch, especially Sec. (4.4) [37]. 
For now, we defer this question for future study [55] . 



B. Lattice Results 

In practice, our method is not quite so trivial. The 

difficulty is that if wc require the comparison lattices to 
have the same physical temperature, but different Nt, 
then the deconfining transition occurs at different val- 
ues of the lattice coupling constant. This significantly 
complicates the analysis. 

In the simulations, the Wilson lattice action was used, 
with lattice coupling constant /3 = 6/g^- The number of 
time steps taken were A^t = 4, 6, 8, and 10. The number 
of steps in the spatial direction, Ns, was always kept fixed 
at Ns = 3Nt; we did not study what happens as this ratio 
is varied. The value of the coupling constant at which 
the deconfining transition occurs, /J^, was determined by 
monitoring the peak in the susceptibility of the triplet 
loop, to give the values in Table I. 

By using non-perturbative renormalization [28], the 
relationship between /? and the temperature was found 
to be: 

\og^ = 1.7139 (^-^J- 0.8155 CS'-:^^) 

-Id 



+ 0.6667 (/f -/3d) , 



(62) 



where f3 = f3 — 6 and = /^^ — 6. In terms of physical 
temperature, our lattices varied from « .ST^ to 3Td. 



TABLE I: The lattice coupling constant for the deconfining 
transition, f3d, at different time steps. 



Nt 


0d 


4 


5.690(5) 


6 


5.89(1) 


8 


6.055(6) 


10 


6.201(5) 



The lattice calculation was done using the over relaxed 
Cabibo-Marinari pseudo-heatbath algorithm. Each up- 
date step contained 4 heatbath updates and six ovcrre- 
laxation steps. A measurement was performed every 10 
update steps. In Table II we summarize our statistics in 
each case. Our lattice data for the bare Polyakov loops 
are presented in Tables HI, IV, V, and VI. They are also 
plotted in Figures 2, 3, 4, and 5. 



TABLE II: . 


Nt 


Measurements 


4 


10000 


6 


400 


8 


400 


10 


400 



TABLE 111: Bart- Bolyakov loops tor X, = 1. 



/3 {^:^) 



(4) 



(4) 



5.50 0.0104(22) 

5.60 0.01659(9) 

5.65 0.02456(14) 

5.69 0.0854(5) 

5.70 0.1233(4) 
5.75 0.17803(17) 
5.80 0.19956(15) 
5.90 0.22978(14) 
6.00 0.25301(13) 
6.10 0.27307(13) 
6.20 0.29075(13) 
6.30 0.30695(12) 
6.40 0.32174(12) 



0.0037(9) 

0.003636(19) 

0.003696(19) 

0.00588(4) 
0.00834(4) 
0.01503(4) 
0.01925(5) 
0.02657(5) 
0.03340(5) 
0.04010(6) 
0.04660(6) 
0.05319(6) 
0.05964(6) 



0.0020(6) 

0.002522(19) 

0.002670(21) 

0.00726(6) 

0.01172(6) 

0.02205(6) 

0.02791(6) 

0.03772(6) 

0.04649(6) 

0.05492(6) 

0.06298(7) 

0.07103(7) 

0.07880(7) 



0.0024(5) 

0.002139(11) 

0.002126(11) 

0.002202(12) 

0.002254(12) 

0.002372(12) 

0.002490(13) 

0.002825(15) 

0.003328(17) 

0.003979(19) 

0.004747(20) 

0.005748(22) 

0.006821(24) 



As the relationship between the lattice coupling con- 
stant and the physical temperature, (62) is nonlinear, it 
is not automatic ensuring that the temperature is the 
same when Nt changes. Thus we resort to interpolation, 
measuring the loops on a fixed grid, in (3, for each Nt. 
For Nt = 4, 6, and 8 wc have linearly interpolated the 
Polyakov loop values to the T/Ta values at which the 
measurements for Nt = 10 were done. Then, for each 
value of T/Td; the expectation value of the bare Polyakov 
loop was fit to (60) — (61). 
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TABLE IV: Bare Polyakov loops for Nt = 6. 



TABLE VI: Bare Polyakov loops for Nt = 10. 



/3 (is) 



{£e} 



s) 



3> 



1 (4> 



6> 



8> 



5.70 
5.80 
5.82 
5.84 
5.86 

5.88 

5.89 

5.90 
5.92 
5.95 
6.00 
6.10 
6.20 
6.30 
6.40 
6.50 
6.60 
6.70 



,00592(16) 

.00963(26) 

.0111(3) 

,0125(3) 

,0183(5) 

,0314(9) 

,0391(10) 

,0545(9) 

,0702(6) 

,0816(5) 

,0935(5) 

,1128(4) 

,1301(4) 

,1445(4) 

,1581(4) 

,1707(4) 

,1829(4) 

,1954(4) 



0.00197(5) 

0.00195(5) 

0.00191(5) 

0.00195(5) 

0.00210(5) 

0.00205(5) 

0.00195(5) 

0.00217(6) 

0.00254(6) 

0.00268(7) 

0.00346(8) 

0.00480(9) 

0.00654(9) 

0.00837(10) 

0.01031(10) 

0.01241(11) 

0.01462(12) 

0.01718(12) 



0.00131(5) 

0.00143(6) 

0.00139(5) 

0.00135(5) 

0.00132(5) 

0.00150(6) 

0.00183(7) 

0.00231(8) 

0.00296(9) 

0.00397(9) 

0.00508(10) 

0.00759(11) 

0.01038(11) 

0.01314(11) 

0.01604(13) 

0.01901(12) 

0.02216(14) 

0.02565(14) 



0.00115(3) 

0.00116(3) 

0.00117(3) 

0.001126(29) 

0.001148(28) 

0.00116(3) 

0.00121(3) 

0.001144(29) 

0.00119(3) 

0.00118(3) 

0.00117(3) 

0.00126(3) 

0.00120(3) 

0.00127(3) 

0.00124(3) 

0.00128(3) 

0.00136(4) 

0.00144(4) 



TABLE V: Bare Polyakov loops for Nt = 8. 



/3 (^3) 



(4) 



(4) 



5.80 
5.90 

6.00 
6.02 0, 
6.04 
6.06 
6.08 
6.10 
6.15 
6.20 
6.30 
6.40 
6.50 
6.60 
6.70 
6.80 
6.90 



00339(9) 
00421(11) 

00668(18) 

,00858(26) 

0159(5) 

0209(5) 

,0345(4) 

0391(3) 

04762(26) 

05400(25) 

,06541(24) 

07540(24) 

08501(22) 

,09492(24) 

10490(24) 

11377(25) 

12214(24) 



0.00133(4 
0.00122(3 
0.00124(3 
0.00126(4 
0.00121(3 
0.00129(3 
0.00129(3 
0.00126(4; 
0.00138(4 
0.00146(4 
0.00164(5 
0.00194(5 
0.00249(5 
0.00299(5 
0.00375(5 
0.00450(6 
0.00540(5 



0.00083(3 
0.00086(3 
0.00087(3 
0.00082(3 
0.00095(4 
0.00092(3 
0.00099(4 
0.00108(4 
0.00131(5 
0.00164(5 
0.00228(6 
0.00305(6 
0.00407(5 
0.00502(6 
0.00635(6 
0.00757(6 
0.00886(7 



0.000775(20) 
0.000770(20) 
0.000763(20) 
0.000775(20) 
0.000754(19) 
0.000713(17) 
0.000751(20) 
0.000755(20) 
0.000754(20) 
0.000724(20) 
0.000768(21) 
0.000751(19) 
0.000804(18) 
0.000765(19) 
0.000766(21) 
0.000819(20) 
0.000778(20) 



In our measurements, we see no statistically significant 
term ~ /tj* ~ 0. Such terms will presumably be 

revealed by more precise mcasiircmonts. The success of 
a fit to (60) indicates that on the lattice, the divergent 
mass m*" does exponentiate. 

We stress that wc make no assumptions about any of 
the functions fn- At a given value of the temperature, 
the logarithm of the expectation value of the bare loop, 
\og({£n)), is a power scries in 1/Nt, beginning as ~ A^t, 
(60) — (61). Given the lattice data, for loops at the same 



6.00 
6.10 
6.15 
6.18 
6.20 
6.22 
6.25 
6.30 
6.40 
6.50 
6.60 
6.70 
6.80 
6.90 
6.95 



00272(9) 

00362(13 

00524(15 

00672(16 

00839(27 

01651(22 

02167(17 

02741(22 

03444(24 

04142(24 

04861(18; 

05505(18 

06204(19 

06852(19 

07130(18 



0.00088(3) 
0.00085(3) 
0.000904(24) 
0.000906(18) 
0.000909(23) 
0.000927(18) 
0.000929(18) 
0.00091(3) 
0.00092(4) 
0.00092(4) 
0.00098(3) 
0.00118(3) 
0.00128(4) 
0.00154(3) 
0.00152(3) 



00058(3) 0. 

00064(4) 0. 

000588(23) 0. 

000615(17) 0. 

000633(24) 0. 

000627(18) 0. 

000641(18) 0. 

00070(4) 0. 

00084(4) 0. 

00099(4) 0. 

00111(4) 0. 

00146(5) 0. 

00193(5) 0. 

00255(5) 0. 

00270(4) 0. 



000525(18) 
000537(22) 
000560(13) 
000550(11) 
000512(14) 
000561(11) 
000532(11) 
000537(18) 
000555(21) 
000521(20) 
000543(17) 
000511(17) 
000549(16) 
000555(16) 
000546(15) 



physical temperature and different values of A^f, there is 
nothing left over to adjust. 

In Fig. 6 wc present a typical fit for the expectation 
value of the bare triplet loop. It is clear that the bare 
loop decreases, with increasing N^, for all temperatures 
measured. 

In Table VII, we give the expectation values of the bare 
Polyakov loops for different representations. Wc chose 
the smallest Nt, Nt = 4, where the signals are greatest. 
For reference, we also include the Casimirs of the different 
representations. 



TABLE VII: Approximate expectation values of the bare 
Polyakov loop, and Casimirs, for Nt = 4 and T = 2Td- 



n 


{in) 




3 


.25 


4/3 


8 


.04 


3 


6 


.035 


10/3 


10 


.004 


6 



For all representations, the signal decreases with in- 
creasing A^^. This indicates that every am*" is positive, 
as suggested before, (57). We were only able to mea- 
sure a signal for the decuplet loop on the smallest lattice, 
Nt — 4. Perhaps a modified Polyakov loop, as discussed 
in the Appendix, might help. We do not discuss the de- 
cuplet loop further. 

Fig. 7 shows our results for the product of the lat- 
tice spacing times the divergent mass, am^". For the 
triplet loop, this product doesn't vary much with tem- 
perature, but those for the sextet and octet loops do. At 
the highest temperature, ~ ST^, the values are approx- 
imately what one expects from lowest order in pertur- 
bation theory, where the divergent masses scale like the 



Casimir of the representation, arrt 



div 

n 



C'n, (50). As 
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the temperature decreases below « 1.5Td, though, the 
perturbative ordering of am*" > am'^™ is reversed, with 
amf''" > am.Q^'" . All divergent masses are approximately 
equal below Td, although the signals are poor. 

After dividing by the renormalization constant, we ob- 
tain the renormalized loops of Fig. 8. Below T^, the 
triplet and sextet fields, which carry Z{?>) charge, should 
vanish in the infinite volume limit. What is striking is 
that below T^, the Z{i) neutral adjoint loop — which 
could be nonzero — is also too small for us to measure: 

(?8)~0 , T<Td. (63) 

Previous studies found this for bare adjoint loops on 
small lattices, such as A^* = 4 [10, 27, 54, 56]. 

In the deconfined phase, the triplet loop is always 
greatest, followed by the octet, and then the sextet loop. 
At T^, the triplet loop jumps to a relatively large value, 
w .4, although the exact value is not very well deter- 
mined. Due to the increase in correlation lengths, there 
is critical slowing down near T^, and much more careful 
studies are required. 

We then computed the expectation value of the differ- 
ence between the sextet and octet loops, and their large 
N limit, (8) and (9). The results are presented in Fig. 
(9). Numerically, we find that both difference loops are 
negative. They each look like a "spike" down, with a 
maximum near T^. The sextet spike is larger, \5^Q\ < .25, 
with a maximum at « 1.25Td. The octet spike is smaller, 
l^^sl l£ -2, with a maximum nearer Td, at w I-IT^. The 
octet spike falls off much more quickly than the sextet 
spike with increasing temperature: by w 1.5Td, the octet 
difference vanishes, while the sextet difference persists all 
of the way to « ST^. 

The most important feature of the difference loops is 
their overall magnitude: each is significantly smaller than 
one. As discussed in the Introduction, this indicates that 
factorization, which is exact for Af = oo, is approximately 
correct for A^ = 3. The difference loops will be discussed 
further in the next section. 

We conclude this section by contrasting our method 
for determining renormalized Polyakov loops with that 
of Kaczmarek et al. [39]. They measure (%(a;)^3(0)) — 
[(£3)1^, and extract by comparing with perturbation 
theory at short distances. Thus they need to measure 
two-point functions, although at just one value of Ni. Wc 
only need to measure one-point functions, but must do so 
at several values of Nf. Our results agree approximately 
with theirs; near T^, we both have |(^3)| « .4. It disagrees 
in that by 3rd, their [(fg)! r:! 1.0, while ours is « .9. We 
do not know the reason for the differences between our 
results, although it could simply be due to the effects of 
finite lattice spacing. 



IV. MATRIX MODELS FOR RENORMALIZED 
LOOPS 

A. Effective Theories 

We now discuss various effective models for the decon- 

fining phase transition [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 
16, 17, 18, 19, 20, 21, 22, 23, 43, 44, 45, 46, 47, 48, 49, 50, 
51, 54], and fit the renormalized triplet loop to a simple 
matrix model [47, 49, 54]. We then discuss the Gross- 
Witten point at large A^, and how it may be related to 
three colors. 

Numerical simulations on the lattice suggest that the 
transition is of second order for two colors [56] and of 
first order for three [24], four [26, 27], and six [27] colors. 
Lucini, Teper, and Wenger [27] present evidence that the 
latent heat grows ~ N"^ for these values of N. From this 
we presume that the transition is of first order for all 
N >2i. Arguments for a first order transition at infinite 
N have also been given by Gocksch and Neri [42]. 

We start with an effective lattice theory in the purely 
spatial dimensions. The actual value of the lattice spac- 
ing is irrelevant: all that matters is that it is much smaller 
than any physical length scale. Thus we concentrate on 
the region near Td, and ignore T — s- 00. 

The simplest approach is to follow Svetitsky and Yaffe 
[5], and construct an effective theory just for the Polyakov 
loop in the fundamental representation. This is certainly 
fine if the transition is of second order, as the only crit- 
ical field should be that for the fundamental loop. For 
the second order transition of two colors, this predicts 
that the universality class is that of the Ising model, in 
agreement with the observed critical exponents [56] . For 
three colors, because Z{2i) symmetry allows a cubic in- 
variant, this approach also predicts that the transition 
is invariably of first order. For more than three colors, 
though, it is not clear why the transition should be of 
first order. It would be if the interactions of the Z[N) 
spins were like those of a Potts model, but they are not at 
high temperature [21]. Loops in higher representations 
can be introduced into this model, as new types of Z[N) 
spins. However, at large N factorization does not follow 
naturally, but must be imposed by hand. 

A variant of this approach is the Polyakov loop model 
[18, 19, 20, 21]. This starts with a potential for the 
(renormalized) fundamental loop, and assumes that the 
pressure is correlated to the value of the potential at its 
minimum. We presume that the same holds for the po- 
tential of the matrix model, sees. IV. B and IV. C. 

The closest model to the underlying gauge theory is a 
gauged, nonlinear sigma model for the Wilson line Ljv 
[6, 18]. In such a model, Higgs phases, in which the 
SU{N) symmetry is spontaneously broken, can appear, 
but these are not expected to arise if there arc no dy- 
namical scalars about. Remember that deconfinement is 
not a Higgs effect: only Z{N) is broken in the deconfined 
phase, not SU{N). 

Assuming that Wilson lines form the essential degrees 
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of freedom in the effective theory, we take as the partition 
function: 



(64) 



i denotes lattice sites in the spatial directions. We assume 
that only Polyakov loops enter into the effective action, 
and consider only constant solutions in mean field the- 
ory. We neglect all kinetic terms, including those for the 
Wilson line, Polyakov loops, and color magnetic fields 
[17]. The effects of fluctuations, which are controlled by 
these kinetic terms, can be important, especially in three 
spacetime dimensions [60]. 

At each site, the measure in (64) implicitly contains 
a constraint to enforce that LAr(i) is a SU{N) matrix. 
As such, loops i-jiii) arc constructed by the usual rela- 
tions of group theory. This is most convenient, since then 
all loops automatically have the right Z{N) charge, and 
satisfy factorization at large A''. Moreover, at each site 
we can use the character expansion, (59), to reduce any 
product of loops to a linear sum. This vastly restricts 
the number of possible couplings which can arise. 

In a sigma model over a symmetric space, the trace 
of L]v is everywhere some fixed constant [38]. In the 
present instance, however, the trace of the Wilson line 
is not constant. Thus the action also includes a poten- 
tial for the Wilson line [18]. Requiring the action to be 
Z{N) invariant, this is a sum over Polyakov loops with 
vanishing Z{N) charge: 



W = A^ S, 



(65) 



the are coupling constants. This series begins with 
the adjoint loop. 

Next, loops on one site can interact with those on an- 
other site: 



A — 

5k = -— Si,ft E^j^+f^'-" pn,n' Re inmwii + n) . 

(66) 

Re denotes the real part, and I3n,w = Pw,tz- We just 
write nearest neighbor interactions, given by the sum 
over the lattice vector fi, but this is inessential. The sum 
over representations is restricted by the requirement that 
the total Z{N) charge of each term is zero, modulo A. 
There arc both diagonal couplings, where TZ' is the rep- 
resentation conjugate to TZ, and off-diagonal couplings, 
where TZ' ^TZ*. This is the complete set of independent 
couplings. 



B. Mean Field and Matrix Models for A = 3 

For three colors, the simplest action includes just the 
triplet loop [48] 



Ss = -3/33 ^i,n Re h{i)tS + n) ; 



(67) 



/33.3» ^ Pz- We use this to develop a mean field approx- 
imation, replacing all nearest neighbors by an average 



value. On a cubic lattice in three dimensions, if the value 
of each neighboring spin is £q = {£3) in (67), the partition 
function of (64) reduces to one for a single site. 



-I 



dLs exp i+lSPsio^eis) = exp(-9V) . (68) 



This is a matrix model, but one whose coupling constant 
depends on the value of the condensate, £q. We introduce 
the single site potential, V. The mean field condition 
is that the average value, computed with this action, is 
equal to the assumed value [38]: 



^0 = -7^ . 
2/33 d£o 



(69) 



As Ljv is an SU{N) matrix, it is the unitary trans- 
formation of a diagonal matrix, with A — 1 independent 
eigenvalues. When A = 2, the integral like (68) is ele- 
mentary, and can be evaluated in terms of Bessel func- 
tions. These can also be done when A > 2, but we found 
it easier to simply evaluate it numerically. Explicitly, 
with L3 = U diag(exp(i0i), exp(j(?2), cxp(— 1(6*1 -|-6'2)) ?7^, 
the normalized Haar measure, including the van der 
Monde determinant, is 



cos(6'i - 6»2))(1 - cos(26'i -t- 6*2)) 



(1 - cos(6'i + 262)) ddide^ . 



(70) 



This mean field theory was studied in the context of 
the deconfinement transition by several groups [47, 48, 
49, 54]. Like the lattice data, the transition is of sec- 
ond order for A = 2, and first order when A > 3. For 
three colors, Damgaard [54] used eq. (69) to compute the 
expectation value of the triplet loop; expectation values 
for the sextet, adjoint and decuplet loops were then com- 
puted from that. Damgaard compared the results of this 
mean field theory to lattice data for bare Polyakov loops, 
with At = 3, by Markum, Faber, and Meinhart [37], 
finding qualitative agreement. 

We stress that the approximate agreement between 
this mean field theory, and lattice data at small A(, is 
in some sense fortuitous. For small Ai, the Z-ji are not 
much different from one, and so the bare values are not 
far from the rcnormalized values; even so, they are not 
identical. To see this in another way, we computed the 
ratio of the difference loop, to the loop itself, for the bare 
octet loop: |((4) - |(4)P)/(4)I. This ratio is - 50% at 
At = 4, and increases to ~ 100 - 200% for At = 10. This 
is to be compared with the values for the renormalized 
octet difference loop in fig. (4), which is < 12%. Thus 
while renormalized loops satisfy factorization, bare loops 
do not. 

In this vein, recently Dittmann, Heinzl, and Wipf com- 
puted the effective potential for bare doublet loops in a 
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pure SU{2) gauge theory [23]. Because only renormal- 
izcd loops satisfy factorization, we suggest that the ef- 
fective potential for renormalized loops is much simpler 
than that for bare loops. 

We next compare the solution of mean field theory, 
eq. (69), to our lattice data for the renormaHzed triplet 
loop. We find that a linear relationship between the mean 
field coupling constant, jSs, and the temperature, T/Td, 
is approximately valid. A least squares fit gives 

Pa = (0.46 ± 0.02) + (0.33 ± 0.02) ^ . (71) 

The computed coupling, as well as the fitted curve, is 
shown in Fig. 10. A quadratic term was included in the 
fit but the coefficient was found to be zero within the 
error bars. Significant deviations are seen at both the 
highest temperature, « BT^, and also at the two points 
closest to Td; see the discussion at the end of this sub- 
section. The approximate linear relationship between (3^ 
and the temperature is typical of mean field theory for 
spin models [38]. 

Using this relationship between the mean field cou- 
pling and the temperature, we then computed mean field 
results for the sextet, octet, and decuplct loops. The 
comparison to our lattice data are shown in Fig. 11. No- 
tice that although we could not extract from the lattice 
a signal for the decuplet loop, the mean field theory pre- 
dicts that while the decuplet loop is less than the sextet, 
it is not that small; for example, ^io(3T<j) Ri .4. 

To obtain a more precise measure of the quality of our 
fits, we computed the difference loops in our mean field 
approximation, and plot them in Fig. 12. We do this 
because even in this simple mean field approximation, 
there are corrections to A/' = oo factorization at A = 3. 
Now compare the difference loops in mean field theory, 
Fig. 12, to those from the lattice. Fig. 9. As expected 
from general arguments, the adjoint difference loop is al- 
ways smaller than the sextet difference loop; also, both 
difference loops are negative in mean field theory. In de- 
tail, however, the difference loops found from mean field 
theory are very different from those found from the lat- 
tice. First, in magnitude the difference loops from mean 
field theory are at least a factor of three times smaller 
than found from the lattice. Further, their tempera- 
ture dependence is very different: in mean field theory, 
both difference loops are greatest at about l.ST^, with 
approximately the same width, ±.5Td. In contrast, the 
difference loops from the lattice have a maximum much 
closer to T^; while the sextet has a tail which persists 
to « 3Td, the octet really appears to be a sharp, narrow 
spike. 

The quality of the fit could be improved by including 
other terms in the action. We started by including an 

adjoint loop in the potential at each site, (65), ^ 78^8- 
Within our numerical accuracy, this only appeared to 
produce a shift in /3s — > /3s + jg. We show in the next 
section that this can be understood at infinite N. 

To model the sextet and octet loops, it is necessary 



to add corresponding fields at each site. In (66) there 
are two diagonal couplings, /^e.e* and P^.h- '^nd one off- 
diagonal coupling, Psfi. Even in mean field theory, it is 
tedious to solve numerically for several, coupled conden- 
sates; the results of a more careful study will be presented 
separately [55]. 

The presence of other loops obviously feeds back into 
the triplet loop. One of the clearest tests of this is the 
value of the triplet loop at T^. In a mean field theory 
which includes just the triplet loop, (69), numerically we 
estimate that is{Td) ~ .485± .001; this is consistent with 
the value of « .49 in [47]. This is significantly higher 
than the values of the renormalized loop from the lattice, 
where both we and [39] find IsiTd) « .4. In a A = 3 ma- 
trix model which includes both triplet and sextet loops, 
we find that the value of the triplet loop at Td decreases 
significantly by adding Ps,6, as this represents a linear 
coupling between the two loops [55] . It is also possible to 
decrease isiTd) by adding a term for the decuplet loop 
to the potential W with the appropriate sign [55]. 

Thus the approximate linearity in /3s with T, (71), 
should be treated as preliminary. Further, we doubt that 
it is true for all couplings. In particular, since the (K;tet 
difference loop is such a sharp, narrow spike in tempera- 
ture, it appears that we can only model it with a /?8 which 
varies non-monotonically with temperature; i.e., which is 
itself a spike. This is not so obvious for the sextet loop, 
due to the triplet-sextet mixing from the coupling Ps,6- 
The sextet loop is also affected by its coupling with itself, 
through the coupling Pe.e, and by the decuplet loop in 
the potential. 

Nevertheless, the matrix model proposed in sec. IV.A 
appears to be a useful way of characterizing the conden- 
sates of renormalized Polyakov loops when A = 3 [55] . 



C. Matrix Models: A > 3 

For general A, the simplest possibility is to start with 
an action including just the fundamental loop, 

A2 

5iv = — ^/3Si,ftRe^jv(i)^^(i + n) , (72) 

/3 = Pn,N'- For the action to be of order A^ at large 
A, /? must be of order one as A ^ oo. Positive val- 
ues of P correspond to a ferromagnetic coupling. As the 
perturbative vacuum at high temperature is completely 
ordered, /3 +oo as T ^ oo. We assume /? > [61]. 

At infinite A, if the fundamental loop condenses, loops 
in higher representations are fixed by factorization. The 
condensate for the fundamental loop is determined by 
the potential, W of (65). 

We start with the mean field analysis of (72), following 
the discussion of Kogut, Snow, and Stone [47]. We then 
discuss how these results change when the potential W, 
(65), is added [45, 46]. A similar discussion was given 
recently by Aharony et al [14]. 
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Replacing the values of all nearest neighbors by an av- 
erage value £o, we need to evaluate the integral at one 
site, 

2 = [ dLjv exp(iV2(2/?4) Re£jv) 



= e^p{-N^VGw{mo)) . (73) 
The mean field condition is 

This condition is equivalent to minimizing the mean field 
potential 

VmfiP,i)^pi^ + VGwm- (75) 

We replace £o by £, and interpret the result as a potential 
for £. At a fixed f3, as usual the vacuum £o is given by 
minimizing Vmf with respect to £. 

Since an overall factor of iV^ is scaled out of the 
potential, a nonzero value of Vmf{£o) implies that the 
free energy ^ A'^^. This is expected in the deconfined 
phase from the liberation of ~ iV^ gluons. In the con- 
fined phase. Thorn [10] noted that as all states are color 
singlets, their free energy is at most of order one. so 
Vm/(^o) ^ ~ 0. This scaling also motivated the 

Polyakov loop model [18, 19, 20, 21]. Here, we assume 
that to go from the single site model to thermodynamics, 
we multiply Vmf by times the volume of space. 

The potential Vgw has been computed in the large N 
limit by Gross and Witten [43]. At infinite N the result 
is nonanalytic, and is given by two different potentials. 
For small £, the potential is just a mass term, 

V-f = f3il-P)£^ . ^<^, (76) 

while at large £, the potential is 

V+f = -2f3£ + f3£' + l log(2/3£) + ^ , £>^. (77) 

The physical interpretation of this potential is rather 
diff'erent from the context in which it arose. Gross and 
Witten considered a U{N) lattice gauge theory in two di- 
mensions, with lattice coupling constant /3gw = 0^ [43]. 
This is the only parameter in the model, and there is 
no condition to fix £o. Instead, the expectation value 
of Re £n is related to the string tension; it changes 
with Pgw, but is always nonzero. The two potentials 
in (76) and (77) correspond to weak and strong coupling 
branches of the free energy. About /3gw = 1/2, the first 
and second derivatives of the free energy arc continuous, 
but the third derivative is not, so there is a third order 
phase transition in Pgw- 

In mean field theory, /3 is an effective coupling for the 
fundamental loop. As a function of £ at fixed /?, the 



first and second derivatives of the potential arc every- 
where continuous, but third (and higher) derivatives are 
discontinuous at a single point, when £ = 1/(2/3). This 
nonanalyticity is special to = cxd: the mean field po- 
tential is everywhere continuous for finite N. 

Overlooking this discontinuity, the potential Vmf be- 
haves as a potential should. When /3 < 1, the potential 
just increases monotonically with £, so the minimum is 
at £q = 0. For (3 > 1, the potential about the origin is 
given by V~y^, and so starts out with a negative mass 
term. The potential decreases with increasing £, with a 
single minimum when 

When £ > £q, the potential increases monotonically. For 
all /3, the potential is bounded at large £, V^j +/3£^ 
as £ oo. (For £ > 1/(2/3). the extremal condition 
dV^f/d£ = is a quadratic equation. There is another 
root besides (78), but it occurs for £ < 1/(2/3), and so 
doesn't matter.) 

Thus there is a confined phase for /3 < 1, and a decon- 
fined phase for /3 > 1. The expectation value of the loop, 
£f), jumps from zero below the transition, /3 = 1~, to i 
just above, (3 = 1+. The latter is (10) of the Introduc- 
tion. In the limit /3 — > oo, Iq ^ 1- 

To verify that the transition is in fact of first order, 
consider the value of the potential at its minimum. In 
the confined phase, £o = 0, so V~y^(^o) = for all /3 < 1, 
including [3 1~ . In the deconfined phase, using (78) 
one finds that 

Kfi^o)^-^ + .-. , (79) 

as /3 1+. Thus the first derivative of Vmf{£o), which 
respect to /?, is discontinuous when (3 = 1. 

If we assume that (3 is linear in the temperature — as 
found for three colors — then the deconfining transition 
is thermodynamically of first order at = oo, with a 
latent heat ~ N^. In this we agree with [47]. 

Even so, the non-analyticity of the potential still has 
striking physical consequences. In particular, exactly at 
/3 = 1, the potential is completely flat for i between 
and 5, = [63]. The potential then increases mono- 
tonically, starting out to cubic order in £ — ^. The only 
reason the order parameter can jump, despite the fiat- 
ness of the potential, is because /3 = 1 is special: then, 
and only then, docs the point at which the potential is 
discontinuous coincide with the nontrivial minimum. 

To appreciate this in another way, consider how the 
mass squared, rri^ = d'^Vmf/d£'^ changes. Approaching 
the transition in the confined phase implies that we com- 
pute about ^0 = 0, 

mi « 2(1 - /3) , /3 ^ 1^ . (80) 
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In contrast, approaching the transition in the deconfined 
phase, we compute about io = ^, 

ml « 4V/? - 1 , (81) 

Thus while both masses vanish at /3 = 1, they vanish 
with different powers of |/3 — 1|. 

The mass of the Polyakov loop is of physical signifi- 
cance in the underlying gauge theory [19, 20]. In coordi- 
nate space, the connected two point function of £{x) is: 

\x\ 

(82) 

In the confined phase, M — <t/T, where a is the 
string tension. In the deconfined phase, one can define 
M = 2mDebye, where moebye is a (gauge- invariant) De- 
bye mass. Assuming that M ^ m, with mass dimen- 
sions made up by some other physical mass scale, such 
as the temperature, this mean field theory predicts that 
the string tension vanishes at the transition as [62]: 

a(T)~(Trf-T)V2 ^ T^T^, (83) 

and the Debye mass, as 

mDebveiT)r^{T-Tay/^ , T^T+. (84) 

One might refer to this as a "critical" first order transi- 
tion: at the transition, the order parameter jumps, but 

the masses vanish, asymmetrically. 

We next include the effects of the potential, (65). We 
start with the simplest term, 72 7^ 0, which is the con- 
tribution of the adjoint loop to the potential. At large 
N, in mean field approximation we need to evaluate the 
integral 

Z= [ dhN exp{N^2f3£Re£N+72\iN\^)) 



= exp(-7V2V(^^,72)) . (85) 
The solution is [45, 46] : 

V{P£, 72) = 72fc' + VawiPe + 12k) ; (86) 

where A; is a variable which one minimizes with respect 
to. The variation with respect to k enforces the condition 

that the expectation value of i'^ satisfies factorization. 
Finally, the mean field solution is given by minimizing a 
potential with respect to t and k: 

Vrnfii, k) - (3e + 72fc2 + VawifM + l2k). (87) 

This is trivial to solve. Expanding about k = £ + Sk, for 
6k = 0, this reduces to the previous mean field potential, 

except that the coupling constant is shifted, /3 ^ /3 + 
72. Further, if £ is extremal with respect to this shifted 



mean field, the term linear in Sk also vanishes. This 
follows because the Gross Witten potential is a function 
only of P£, and not of /3 and £ separately. As discussed 
previously, for = 3 we discovered numerically that in 
mean field theory, we can shift the adjoint coupling away, 
/3^/3 + 72. 

This is very different from lattice gauge theories in two 
dimensions, where /3gw and 72 represent independent 
coupling constants [45, 46]. Then f^aw and 72 can be 
varied irrespectively of each other, and one finds that the 
third order transition, in Paw for 72 = 0, can become a 
first order transition in the plane of (3gw and 72 [45, 46]. 

The generalization to 74 is direct. One includes 
a constraint field for |^|^ and then solves the constraint 
at A'' = cx). The mean field coupling (3 is again shifted, 
/3 ^ /? + 72 + #74^0' by an amount which depends 
on the condensate, £0; there are also additional potential 
terms, ~ 74 ^o- Consequently, 74 ^ corresponds to a 
change in the potential. Aharony et al. [14] show that 
for 74 7^ 0, ordinary transitions appear to be generic. In 
particular, first order transitions have finite correlation 
lengths in both phases. 

This remains the case for an arbitrary potential. W. At 
large A^, this is just a sum of powers of the fundamental 
loop: 

W = Ei^m (72m {\£N\^r + Un ((^iv)""^ + C.C.)) . 

(88) 

This is a typical potential for a scalar field £n = ^n('')- 
the adjoint loop acts like a mass term, while the other 
terms are interactions of quartic and higher order, invari- 
ant mider a global symmetry of Z{N). 

We define the Gross-Witten point as the transition 
for a potential which is just the adjoint loop. All other 
interactions are dropped: /? -|- 72 ^ 0, with 74 = 76 = 

■ ■ ■ = — £,2N — . . . 

The deconfining transition for three colors appears to 
be close to the Gross-Witten point of infinite N . For 
example, although the data are very limited [25] , the de- 
crease of the Debye mass near does seem to be signif- 
icantly sharper than for the string tension, as indicated 
by the different "critical exponent" in (84) versus (83). 
Exactly how close A" = 3 is to the Gross Witten point of 
N = 00 can be characterized within matrix models [55]. 
Effects which are important for three colors include the 
contribution of the decuplet loop, which is like a cubic 
interaction for the triplet loop, and the mixing between 
the triplet and sextet loops, P^.q in (66). 

Assuming that three colors is near the Gross-Witten 
point, we can explain why Z{N) neutral fields have small 
expectation values in the confined phase, (63). In the 
confined phase, the potential is purely a mass term, 
as corrections to (76) begin with the baryon vertex, 
- (^w)^ [44, 47, 51]. This vertex induces Z{N) neutral 
expectation values, but as noted by Goldschmidt [44], 
these are of order ^ exp(— A'). 

That N = 3 is close to the Gross- Witten point could 
be an accident of three colors. The lattice will tell us 
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if the dcconfining transition for four or more colors is 
also close to the Gross Witten point. The lattice finds a 
first order transition [26, 27], but the crucial tests arc the 
value of the renormalized, fundamental loop at T^, and 
whether the string tension and the Debye mass decrease 
significantly near Td- 

If this is not found, the most probable scenario is just 
that the transition becomes more strongly first order with 
increasing N . We term the transition strongly first order 
if the value of the renormalized fundamental loop at is 
near unity. Then at the transition, deconfinement is not 
halfway, as it is at the Gross-Witten point, but nearly 
complete. 

For such a strongly first order transition, neither the 

string tension, nor the Debye mass, would need to change 
much about T^. Gocksch and Neri [42] argued that the 
"loop" string tension is constant for T <Td. In a Nambu 
string model, though, a large d expansion shows that 
when the ordinary string tension vanishes at T^, the loop 
string tension is still « 87% of its value at zero temper- 
ature [62]. We are unaware of lattice data on the loop 
string tension. 



V. OUTLOOK 

In this paper we presented a general analysis of the 
renormalization of Polyakov loops, and applied it to mea- 
sure the simplest loops for three colors in four dimensions. 
These results led us to consider effective matrix models 
for the deconfining phase transition. There are clearly 
many avenues for future study. 

For two colors, there will be two regions. About T^, 
there is a critical region, controlled by universality of the 
second order transition [5], and in which factorization 
fails. This may then match onto a mean field region, 
where factorization is approximately valid. 

For three colors, careful measurements of the renormal- 
ized loops, and associated masses, will sharply constrain 
the couplings of the effective matrix model [55]. We note 
that while we could not extract an expectation value for 
the renormalized decuplet loop from that for the bare 
loop, mean field theory indicates that it is significant. 

Simulations should quickly show if for the deconfin- 
ing transition for four or more colors is near the Gross- 
Witten point. 

Considering theories other than SU{N), Holland, 
Minkowski, Pepe, and Wiese [65] noted that in a pure 
G(2) gauge theory, there is no center to the gauge group, 
and so no absolute measure of confinement. This is analo- 
gous, though, to Z{N) neutral loops in SU (N), for which 
we measured no signal below T^. In the simplest mean 
field theory for a G{2) gauge theory, presumably there is 
a first order transition, with a value for the fundamental 
loop at near i. Maybe like SU{3), the "deconfin- 
ing" transition in a G{2) gauge theory is also near the 
Gross-Witten point. 

The renormalization of Wilson lines implies that once 



the divergent mass is known, all renormalized loops can 
be computed. We suggest that numerical simulations 
measure loops of different shapes, such as Polyakov loops 
with cusps. Fig. (1), and circular loops, as can be com- 
puted in supersymmetric theories [13]. 

In the end, however, what is most important is to 
measure renormalized Polyakov loops for theories with 
dynamical quarks. Our method for computing renormal- 
ized Polyakov loops is completely unaffected by the pres- 
ence of dynamical quarks. Given the flavor independence 
found for the pressure [24], it would be striking if the 
values of renormalized Polyakov loops, with dynamical 
quarks, are found to be close to those of the pure gauge 
theory. For recent results, see [39]. 
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APPENDIX: IMPROVED WILSON LINE 

It is difficult extracting renormalized Polyakov loops 

in representations such as the decuplet, because the bare 
loop is suppressed by a small renormalization constant. 
In this Appendix we give a formal discussion of how to 
improve the Wilson line [57] . 

Our discussion applies to any Wilson line along a path 
x^{s), where s is the path length along the curve. As 
shown in sec. II, (13), the propagator for a test quark is 
proportional to the Wilson line. Consequently, we con- 
sider a generalized propagator, by adding an operator X 
to the covariant derivative: 

(^^ - igA'^x^' - g{s, s') = 6{s - s') ; (A.l) 

where x^^ = dx^^/ds. The representation is denoted im- 
plicitly. Schematically, the solution to this equation is 



g{s, s') = 0{s-s') V exp {igAi'xi' + X)dsj , (A.2) 

with V denoting path ordering. The solution to Q is 
schematic because of the path ordering, but it is easy to 
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understand the sohition as a power scries in X, with each 
insertion of X sandwiched between a Wilson line on both 
sides. 

Any possible operator X has a higher mass dimension 
than the gauge field, so we make up the dimensions with 
inverse powers of the ultraviolet cutoff, such as the lat- 
tice spacing a. The important thing is that X respects 
the relevant symmetries. Gauge invariance requires that 
X transforms homogeneously under gauge transforma- 
tions, but that is simply done by using powers of the 
field strength tensor, G^,^. The operator must also be in- 
variant under how we parameterize arc length, s s'(s). 
The final symmetry is the zig-zag symmetry of Polyakov 
[66] ; for reparameterizations which go in the opposite di- 
rection, with ds/ds' < 0, A" should change sign. 

The simplest possibility is 

X = a'g^ G%Vi^ . (A.3) 

This is gauge covariant and reparameterization invariant, 
but is not zig-zag symmetric. If one abandons zig-zag 
symmetry, then this operator can be used to regularize 
the Wilson line, since it is just a field strength dependent 
"mass" term for the line. 

One can satisfy all symmetries with the following op- 
erator, 

X = agG^,^—^ , (A.4) 

where x = d'^x'^ /ds^. A similar operator was noted 
by Polyakov and Rytchkov [66]. Since it vanishes for 
a straight path, where x^ = Q, \t doesn't help to regu- 
larize the Wilson line. It is also problematic to continue 
to Minkowski spacetime, since it is singular on the light 
cone [66]. 

To define an operator which satisfies aU of our require- 
ments, we define a unit vector normal to the path, 

h^x^" = , = 1 . (A.5) 

For a given direction, we introduce the h dependent op- 
erator 

Xn = Kg G^^x'^h" . (A.6) 
We then define a modified Wilson line as 

J dClf, V exp (^J {igA^x^^ + Xf,)d^ . {A.7) 

The operator inserts the field strength tensor perpen- 
dicular to the path. We then integrate over all directions 
of the insertion, with dD,n the normalized integral over 
n, J dilfi = 1. We obviously cannot integrate over all di- 
rections n in the exponential, or the term would vanish, 
and so do so in the prcfactor of the exponential. Zig-zag 
symmetry is maintained by integrating over all n. 

In perturbation theory, to lowest order Xn generates 
new divergences. In three spatial dimensions, and in- 
cluding a Debye mass tod ~ gT in the propagator for Aq 



[3] , the leading divergence is 

We assume that the ultraviolet divergence is cutoff 
strictly at momenta 1 /a. The usual term is ~ / d^k/k^ ~ 
1/a, and is of the same order as the now term, ~ 
a^/ d^k ~ 1/a. For k < 1/a/3, (A.8) is negative, while 
for K > 1/ -\/3, it is positive. If a more realistic cutoff is 
used, then the value of n at which the sign changes will 
be different, but for large k (and to leading order in g"^) 
bare loops are enhanced, not suppressed. 

Following the procedure in sec. (II. D), from log((£)) 
we compute as power scries in Nt- The term linear in Nt 
gives Z-fi, with the renormalized loop given by the term 
independent of Nt- We argue that while the Zti are k 
dependent, the renormalized loops are not, at least in 
perturbation theory. To lowest order in g^, for k = 
the renormalized loop arises from the correction from 
the Debye mass term, ~ g^{m%y/'^/T ~ g^, (55) [3]. 
This is non-analytic in the Debye mass, and arises be- 
cause the leading term is only linearly divergent. For 
the K-dependent term, the leading term is cubically di- 
vergent, so corrections are ~ g'^K?'a?m\,/{aT) ~ g'^n^aT; 
but this is ^ 1/Nt, and vanishes as Ni oo. There is a 
non-analytic term at one higher order in momentum, but 
this is even smaller, ~ g'^ a? {m\,Y/'^ /T ~ g^K^{aT)^ 
l/Nf. It seems likely that this holds for any K-dependent 
terms: they contribute to terms ~ l/{aT) ~ TVj, or to 
lattice corrections ~ 1/Nt or smaller, but not to terms 
in the continuum limit, ~ N^ . This analysis is special 
to four spacetime dimensions: in three dimensions, the k 
dependent term is quadratically divergent, and has loga- 
rithmic corrections, as for k = 0. 

While zig-zag symmetric, the term added, Aft, is not 
a phase factor; the coupling k must be real in order to 
enhance the bare loop. As the bare loop is not the trace of 
a unitary matrix, there is no bound on the renormalized 
loop, as for k = 0, (56); when 0, 2^ can diverge in 
the continuum limit, instead of vanishing. 

On the lattice, adding a field strength tensor at a point 
corresponds to stapling a plaquctte, as the average over n 
becomes a sum over the directions transverse to the path. 
To lowest order in k, this modification is the same as the 
smearing of link variables proposed by the APE collab- 
oration [57]. We suggest doing this not just to lowest 
order, but to all orders in k. This is related to the stout 
links of Morningstar and Peardon [57] . The difficulty is 
that eventually the Wilson line is smeared over the entire 
lattice, which must then be cut off in some way. On the 
other hand, the usual problem with smeared links is the 
need to project back to an element of SU{N), which is 
unnecessary here. 
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FIG. 1: A Polyakov loop with two cusps, at r = and r — 
1/{2T). The dotted hnes denote t = and l/T. 
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FIG. 2: The bare triplet Polyakov loop as a function of tem- 
perature. 
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FIG. 3: The bare sextet Polyakov loop as a function of tem- 
perature. 
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FIG. 4: The bare adjoint Polyakov loop as a function of tem- 
perature. 
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FIG. 5: The bare decuplet Polyakov loop as a function of 
temperature. 
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FIG. 6: The bare triplet Polyakov loop at fixed temperature, 
as a function of Nt- 
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FIG. 7; The divergent masses, times lattice spacing, versus 
temperature. 




FIG. 8: Renormalized Polyakov loops as a function of tem- 
perature. 
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FIG. 9: Difference loops: renormalized Polyakov loops, minus 
their large values. 




FIG. 10: The fundamental coupling constant extracted from 
the lattice data. The circles are the fundamental coupling 
computed from lattice data for the fundamental loop. The er- 
ror bars on the extracted points are smaller than data points. 
The smooth line is the linear fit to the extracted fundamental 
coupling with the errors quoted in (71). 
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FIG. 11: The values of the Polyakov loops from the matrix 
model, using the linear relationship between the coupling and 
the temperature. The lower-most (green) line corresponds to 
the value of the decuplet loop, for which there is no lattice 
data. 




FIG. 12: Difference loops from the matrix model: Polyakov 
loops, minus their large values. 



